Self-dual Lorentzian Wormholes and Energy in Teleparallel Theory of Gravity * 
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Two spherically symmetric, static Lorentzian wormholes are obtained in tetrad theory of grav- 
itation as a solution of the equation p = p t = 0, where p = T^rtu 3 , p t = (Ty — ^Tgij)u % u 3 and 
u l Ui = — 1. This equation characterizes a class of spacetime which are "self-dual" (in the sense of 
electrogravity duality). The obtained solutions are characterized by two-parameters k\, k-i and 
have a common property that they reproduce the same metric spacetime. This metric is the static 
Lorentzian wormhole and it includes the Schwarzschild black hole. Calculating the energy content 
of these tetrad fields using the superpotential method given by M0ller in the context of teleparallel 
spacetime we find that E = m or 2m which does not depend on the two parameters k\ and k^ 
characterize the wormhole. 
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1. Introduction 



It was recognized by Flamm [1] in (1916) that our universe may not be simply connected, there 
may exist handles or tunnels now called wormholes, in the spacetime topology linking widely sepa- 
rated regions of our universe or even connected us with different universes altogether, traversable 
Lorentzian wormholes "which have no horizons, allowing two-way passage through them", and 
were especially stimulated by the pioneering works of Morris, Thorne and Yurtsever [2, 3], where 
static, spherically symmetric Lorentzian wormholes were defined and considered to be an exciting 
possibility for constructing time machine models with these exotic objects, for backward time 
travel [4, 5, 6]. 

Morris and Thorne (MT) wormholes are static and spherically symmetric and connect asymp- 
totically flat spacetimes. The metric of this wormhole is given by 

r\r 2 

ds 2 = _ e 2Hr) df 2 + + r 2( d9 2 + ^2 ^2) (1) 

1 — b(r)/r 

where $(r) being the redshift function and b(r) is the shape function. The shape function describes 
the spatial shape of the wormhole when viewed. The metric (1) is spherically symmetric and 
static. The coordinate r is nonmonotonic in that it decreases from +00 to a minimum value bo, 
representing the location of the throat of the wormhole, and then it increases from 6 to +00. 
This behavior of the radial coordinate reflects the fact that the wormhole connects two separate 
external universes. At the throat r = b = b , there is a coordinate singularity where the metric 
coefficient g rr becomes divergent but the radial proper distance 

l{r) = ±r , ^ , (2) 
•>b yjl-b(r)/r 

must be required to be finite everywhere [7]. At the throat, l(r) = 0, while l(r) < on the 
left side of the throat and l{r) > on the right side. For a wormhole to be traversable it must 
have no horizon which implies that gu must never allowed to be vanish, i.e., <J>(r) must be finite 
everywhere. 

Traversable Lorentzian have been in vogue ever since Morris, Thorn and Yurtsever [3] came up 
with the exciting possibility of constructing time machine models with these exotic objects. (MT) 
paper demonstrated that the matter required to support such spacetimes necessarily violates the 
null energy condition. Semiclassical calculations based on techniques of quantum fields in curved 
spacetime, as well as an old theorem of Epstein et. al. [8], raised hopes about generation of such 
spacetimes through quantum stresses. 

There have been innumerable attempts at solving the "exotic matter problem" in wormhole 
physics in the last few years [5, 9]. Alternative theories of gravity [10] evolving wormhole space- 
times [11]~ [14] with varying definitions of the throat have been tried out as possible avenues of 
resolution. 

M0ller modified general relativity by constructing a new field theory in Weitzenbock spacetime 
[15]. The aim of this theory was to overcome the problem of energy- momentum complex that 
appears in the Riemannian spacetime [16]. The field equations in this new theory were derived 
from a Lagrangian which is not invariant under local tetrad rotation. Sdez [17] generalized M0ller 
theory into a scalar tetrad theory of gravitation. Meyer [18] showed that M0ller theory is a special 



The tetrad theory of gravitation based on the geometry of absolute parallelism [21]~[28] can 
be considered as the closest alternative to general relativity, and it has a number of attractive 
features both from the geometrical and physical viewpoints. Absolute parallelism is naturally 
formulated by gauging spacetime translations and underlain by the Weitzenbock spacetime, which 
is characterized by the metric condition and by the vanishing of the curvature tensor using the 
affine connection resulting from the geometry of the Weitzenbock spacetime. Translations are 
closely related to the group of general coordinate transformations which underlies general relativ- 
ity. Therefore, the energy-momentum tensor represents the matter source in the field equation for 
the gravitational field just like in general relativity. 

It is the aim of the present work to derive a spherically symmetric solutions in the tetrad 
theory of gravitation. To do so we first begin with a tetrad having spherical symmetry with three 
unknown functions of the radial coordinate [29]. Applying this tetrad to the field equations of 
M0ller's theory we obtain a set of non linear partial differential equations. We propose a specific 
restriction on the form of the stress-energy that when solving the non linear partial differential 
equations, automatically leads to a class of wormhole solutions. The characterization of our class 
of self-dual wormholes is 

P = Pt = 0, (3) 

where p and p t are respectively the energy density measured by a static observer and the conver- 
gence density felt by a timelike congruence. Indeed Eq. (3) characterizes a class of "self-dual" 
wormhole spacetimes which contains the Schwarzschild solution. The notation of duality involves 
the interchange of the active and passive electric parts of the Riemann tensor (termed as electro- 
gravity duality) [13]. Electrogravity duality essentially implies the interchange of the Ricci and 
Einstein tensors. For vanishing Ricci scalar tensor the Ricci and Einstein tensors become equal 
and the corresponding solution could be called "self-dual" in this sense. 

The weak energy condition says that the energy density of any system at any point of spacetime 
for any timelike observer is positive (in the frame of the matter this amounts to p > and p+p > 0). 
When the observer moves at the speed of light it has a well defined limit, called the null energy 
condition (in the frame of matter p + p > 0) [30]. 

Under the duality transformation, p and p t are interchanged indicating invariance of Eq. (3). 
Since energy densities vanish and yet the spacetime is not entirely empty, the matter distribution 
would naturally have to be exotic (violating the energy condition) [13]. Physically, the existence of 
such spacetimes might be doubted because of this violation of the weak and null energy conditions. 
Analogous to the spatial-Schwarzschild wormhole, for which g 00 = —1 and gu = (1 — , 
these spacetimes have zero energy density but nonzero pressures [13]. The spatial-Schwarzschild 
wormhole is one specific particular solution of the equations p — 0, p t — 0. The central problem of 
the traversable wormhole is connected with the unavoidable violation of the null energy condition. 
This means that the matter which should be a source of this object has to possess some exotic 
properties. For this reason the traversable wormhole cannot be represented as a self-consistence 
solution of Einstein's equations with the usual classical matter as a source because the usual matter 
is sure to satisfy all the energy conditions [31]. Therefore, it is a basic fact for the construction of 
the traversable wormholes that the null energy condition has to be violated [30]. 

It is our aim to discuss the physical properties of the spherically symmetric wormholes obtained. 
In section 2 a brief survey of Miller's tetrad theory of gravitation is presented. The exact solutions 
of the set of non linear partial differential equations are given in section 3. In section 4 the energy 
content of these solutions are given. Discussion and conclusion of the obtained results are given 
in section 5. 



2. M0ller's tetrad theory of gravitation 



In a spacetime with absolute parallelism the parallel vector fields define the nonsymmetric 
affine connection 

r A d = e V U) 
where e* Pi ^ = d v e l ^ The curvature tensor defined by T x ^ u is identically vanishing, however. 

M0ller's constructed a gravitational theory based on this spacetime. In this theory the field 
variables are the 16 tetrad components e^, from which the metric tensor is derived by 

<T = ^Ve/, (5) 

where r/ y ' is the Minkowski metric rjij = diag(+l , — 1 , — 1 , —1). 

We note that, associated with any tetrad field there is a metric field defined uniquely 
by (5), while a given metric g^ v does not determine the tetrad field completely; for any local 
Lorentz transformation of the tetrads leads to a new set of tetrads which also satisfy (5). The 
Lagrangian L is an invariant constructed from 7 MJ , P and g^, where 7 (UI/p is the contorsion tensor 
given by 

IfjLUp f] ^ &i p&jv; pi (6) 

where the semicolon denotes covariant differentiation with respect to Christoffel symbols. The 
most general Lagrangian density invariant under the parity operation is given by the form [15] 

C ^ (cn^ + aa'T'V, + «37^ P 7^) , (7) 



where 



g d = det(^), (8) 



and $ M is the basic vector field defined by 

% d = 7V- (9) 

Here ai,a 2 , and 0:3 are constants determined by M0ller such that the theory coincides with general 
relativity in the weak fields: 

«! = --, a 2 = -, a 3 = -(l-2A), (10) 

K K K 

where k is the Einstein constant and A is a free dimensionless parameter*. The same choice of the 
parameters was also obtained by Hayashi and Nakano [22]. 

M0ller applied the action principle to the Lagrangian density (7) and obtained the field equation 
in the form 

G^ + H^^-kT^ K^ = 0, (11) 
where the Einstein tensor G^ u is the Einstein tensor defined by 

Gp V d = f - iWtt) - \g,M{}), 



where i? M ^({}) and R({}) are the Ricci tensor and Ricci scalar. and K^ v are given by 

,, dcf. . r ... ( l 



Ipapl^ v ^ r 1paix1v' + 'paul ' p? ^ 9pv ( 1 pari P „1 porl^ 



(12) 



and 



lpu P . p 



(13) 



and they are symmetric and skew symmetric tensors, respectively. 

M0ller assumed that the energy-momentum tensor of matter fields is symmetric. In the 
Hayashi-Nakano theory, however, the energy-momentum tensor of spin-1/2 fundamental parti- 
cles has non-vanishing antisymmetric part arising from the effects due to intrinsic spin, and the 
right-hand side of antisymmetric field equation (11) does not vanish when we take into account 
the possible effects of intrinsic spin. 

It can be shown [23] that the tensors, H pu and K^, consist of only those terms which are 
linear or quadratic in the axial- vector part of the torsion tensor, a M , defined by 



def. £ 
a p g e pupal 



vpa 



where 



dcf 



^pupa 



lei . I r 



pupcn 



(14) 



where 5 pupa being completely antisymmetric and normalized as #0123 = — 1. Therefore, both H pu 
and vanish if the a p is vanishing. In other words, when the a p is found to vanish from the 
antisymmetric part of the field equations, (11), the symmetric part of Eq. (11) coincides with the 
Einstein field equation in teleparallel equivalent of general relativity. 



3. Spherically Symmetric Solutions 



Let us begin with the tetrad [29] (lines are labelled by I and columns by /i) 



(A 



Dr 



B 



B cos 9 



B 

sin 9 

r 



\ 

B sin d 



B sin 9 cos <\> — cos 9 cos d> 

r r sin 9 

B Bcosd> 

5 sin sin — cos sin — 

r r sin 9 



) 



(15) 



where A, D, B, are functions of the redial coordinate r. The associated metric of the tetrad (14) 
has the form 



ds* = - B \ 2 °T de - 2^drdt + ^ + "-^ + sin 2 9d^). 



(16) 



As is clear from (16) that there is a cross term which can be eliminated by performing the 
coordinate transformation [32] 



im 1 1 



ADr 



using the transformation (17) in the tetrad (15) we obtain 



A 



1 - V 2 R 2 
AVR sin 9 cos 

1 - V 2 R 2 
AVR sin 9 sin 

1 - V 2 R 2 

AVR cos 9 
V 1 - V 2 R 2 



(RV - R 2 VB') 
(1-RB') sin 9 cos 
(1-RB') sin 9 sin 

(1 - RB') cos 9 





cos 9 cos 4> sin 

# "^snTfl 

cos 9 sin cos 



R 

sin # 



i?sm# 




(18) 



where A, V and B are now unknown functions of the new radial coordinates R which is connected 
to the old redial coordinate r through the transformation defined by 

* dB(R) 



R 



B' 



with B' 



dR 



(19) 



We are interested in finding some special solutions to the partial non linear differential equa- 
tions resulting from applying the tetrad (18) to the field equations (11). 

The First Solution 

If the unknown function V(TZ) = 0, then the resulting partial non linear differential equations 
have the form 



p(R) = -^-{3B' 2 R-2RB"(R) + 2R 2 B'(R)B"(R)-4B'(R)}, 
r(R) 



kRA 
p(R) = kT\ 



l —\2B'A- ARA'B' + 2 A' + 2R 2 A'B' 2 - RAB' 2 } , 

[2R 2 A ,2 B' - A[R 2 A'B']' + A 2 [RB']' + A[RA']' - 2RA") , (20) 



where 



p{R) = T%, 



r(R)=T\, 



p{R) = T 2 2 = T 3 



with p(R) being the energy density, r(R) is the radial pressure and p(R) is the tangential pressure. 
(Note that t(R) as defined above is simply the radial pressure p r , and differs by a minus sign from 
the conventions in [2, 5].) 

Eqs. (20) can be solved to take the form 



A(K) 



h + k 2 \ 1 - 



2m 
~R 



B(K)=ln{RlR- m + RJl- 2 ^\\-2Jl- 2 -^, (21) 



where k\ and k 2 are constants of integration. The associated Riemannian metric of solution (21) 
takes the form 



ds 2 = -rn(R)dT 2 ^ 



dR 2 



+R 2 dtt 2 , where Vl (R) = ^ + fc 2 ^l-^j , V2 (R) = (l - ^) , 

(22) 



The Second Solution 

If the unknown function B(1Z) = 1, then the resulting partial non linear differential equations 
have the form 

p(R) = —{3V + 2RV}, 
r(R) = - V{1 " RV ' 2) {3V + 2RV'}, 



A 2 K 



p{R) = kT\ 



— — { R 3 A 2 VV" + 6R 2 A 2 W + 3RA 2 V 2 + AA' - 4R 2 AV 2 A' 
kRA 2 \ 

{2RA r2 - RAA"){1 - R 2 V 2 ) - 3R 3 AVA'V + R 3 A 2 V ,2 \. (23) 



Eqs. (23) can be solved to take the form 



A{K) 



k 2 + 



V(K) 



1 2m 



2m 



(24) 



Using (21) or (24) in (11) we can get the components of the energy-momentum tensor turn out 
to have the form 



p(R) = 0, t(R) 



The weak energy 



and null energy conditions 



2mk\ 



R 3 Ui + hJl 



2m 
~R 



V (R) = - 



mki 



r 3 fci + yi 



2m 



p>0, p + r>0, p + p>0, 



p + r>0, p + p>0 



(25) 
(26) 

(27) 



are both violated as is clear from (25). The violation of the energy condition stems from the 
violation of the inequality p + r > 0. 

The associated Riemannian space of solution (24) has the form (22). If one replacing k 2 by 
— k 2 at the above solutions, ((21) or (24)), the resulting form will also be a solution to the non 
linear partial differential equations ((20) and (23)). The Ricci scalar tensor vanishing, i.e., 

*({})= 0, 

for solutions (21) and (24). The metric (22) makes sense only for R > 2m so to really make the 
wormhole explicit one needs two conditions patches 

i?i G (2m, oo), i?2 £ (2m, oo), 

which we then have to sew together at R = 2m. 



Transforming the metric (22) to the isotropic coordinate using the transformation 



(28) 



Using (28) in (22), then the transformed metric will have the form 



ds 2 = -rj^dT 2 + m (R) 



dR 2 + R 2 dn 2 



(29) 



where r)i(R) and r) 2 (R) have the form 



rn{R) = \k 1 + k 2 



1 



m 
2R 



m 

1 + m 



(30) 



and the energy momentum tensor given by Eq. (25) takes the form 



r(R) = 



128A;im 



R 2 (l + [2R(h + k 2 ) + m(ki - k 2 ) 



p{R) 



QAkim 



K 



2R(ki + k 2 ) + m(h - k 2 ) 



(31) 



As is clear from Eqs. (31) that the energy momentum tensor t(R) and p(R) have a common 
singularity, i.e., R ~ 0. This is due to the advantage of the isotropic coordinates. When the 
geometry is such that it can be interpreted as a Lorentzian wormhole then the isotropic coordinate 
patch is a global coordinate patch. Now one has a single coordinate patch for the traverse wormhole 
which will be use to discuss the properties of the geometry. 

Let us explain this result by studying the geometry of each solution. 
1) The geometry of the two solutions (21) and (24) is invariant under simultaneous sign flip 



k 2 — > ~k 2 , 



it is also invariant under simultaneous inversion 

o 

tS m i 

R — > — ^ and sign reversal 

keeping k\ unchanged. 
2) 

k 2 ^0, k ± = 0, 

gives the Schwarzschild geometry, it is non-traversable. 

3) 



k 2 



-k 2 , 



k 2 = 0, 



4) 

k 2 = 0, fci = 

is a singular. 

5) At the throat g tt (r = 2m) = — fc 2 , so fci ^ is required to ensure traversability. 

6) To see if there is ever horizon or not Dadhich et. al. [13] done the following discussion. 

A horizon would seem to form if the component g tt = 0, i.e., if there is a physically valid 
solution to 

*i(l + |g)+*2(l-|g)=0, (32) 
- _ m k 2 - fci 

* fc " TfcTTfc? (33) 

* 2 - ^ > 0, (34) 



whose solution has the form 



that is a horizon tries to form if 



k 2 + h 

which is actually a naked singularity. This singularity occurs if either 

k 2 + ki > 0, and k 2 — k\ > 0, 

or 

k 2 + ki < 0, and k 2 — k\ < 0. (35) 

Outside of these regions the naked singularity does not form and one has a traversable wormhole. 
Dadhich et. al. [13] studied the k\ — k 2 plan and have found that if 

k 2 = Z cos 9, ki — Z sin 9, 

i) 9 = 0: Schwarzschild spacetime. 

ii) 9 G (0, 7r/4), naked singularity. 

hi) 9 = 7r/4 => k 2 = ki, i.e., the R — > region is not flat. 

iv) 9 G (7r/4, 37r/4) traversable wormhole. 

v) 9 = 7r/2 is the spatial-Schwarzschild wormhole. 

vi) 9 = 3n/A =^ k 2 = —fci, the R — > oo region is not flat. 

vii) 9 G (37r/4, 7r) naked singularity. 

viii) 9 = 7r: Schwarzschild spacetime. 

ix) 9 > n: repeat the previous treatment. 

The ratio value of the two parameters fci and k 2 that characterize the solutions (21), (24) is 



obtained fci/fc 2 > 10 8 7 with 7 2 /5 2 < 2 implying (5 < ^2/3 [13] . 

Thus we have two exact solutions of the field equations (11), each of which leads to the same 
metric as given by Eq. (22). The axial vector part of the two solutions (21) and (24) vanishing 
identically i.e., a M = 0. Therefore, the two tensors H^ u and K^ u are also vanishing and Miller's 
tetrad theory coincides in that case with the teleparallel equivalent of general relativity. 



4. Energy content 



The superpotential is given by 

U P vX = i - Z £^ P xP™ X l<f> P 9 ax 9,r - X9r,l xpa - (1 - 2X)g T ,r px ] , (36) 

where P xp<J TvX is 

Pxp« TuX =■ 8x9p* vX + 5 P T 9«x X ~ 8* T 9 XP vX (37) 
with g pa vX being a tensor defined by 

9 P a vX d = 5 p »5 a x - 5/5 p \ (38) 



The energy-momentum density T p is defined by 



j, dcf . i , u\ 
M — M M , A' 

and automatically satisfies the conservation law, r p v = 0. The energy is expressed by the surface 
integral [33, 34, 35] 

E = lim / U 0a n a dS, (39) 

r ^°° j r=constant 

where n a is the unit 3-vector normal to the surface element dS. 

Now we are in a position to calculate the energy associated with solution (21) using the 
superpotential (36). As is clear from (39), the only components which contributes to the energy is 
U 0a . Thus substituting from solution (21) into (36) we obtain the following non-vanishing value 

M °°° = 75T (40) 

Substituting from (40) into (39) we get 

E = m. (41) 

Repeat the same calculations for solution (24) we obtain the necessary components of the super- 
potential 

U ° (42) 

and the energy will have the form 

E = 2m. (43) 



5. Discussion and conculusion 



In this paper we have applied the tetrad having spherical symmetry with three unknown 
functions of radial coordinate [29] to the field equations of M0ller's tetrad theory of gravitation 
[15]. From the resulting partial differential equations we have obtained two exact non vacuum 
solutions. The solutions in general are characterize by three parameters m, k\ and k 2 . If the two 
parameters k\ — and k 2 — \ then one can obtains the previous solutions when the exponential 
term equal zero, i.e., e _ij3//rl3 = [36]. The energy-momentum tensor has the property that 
p(R) = for the two solutions. This leads to the violation of the weak energy and the null energy 
conditions defined by Eqs. (26) and (27) due to the fact that the radial pressure is negative, i.e., 
t(R) as is clear from Eq. (25). The line element associated with these solutions has the same 
form (22). 

To make the picture more clear we discuss the geometry of each solution. The line element of 
these solutions in the isotropic form are given by (29). If gu = one obtains a real naked singularity 
region. Outside these regions naked singularity does not form and one obtains a traverse wormhole. 
The throat of this wormhole gu(R = 2m) gives the conditions that g tt = —k\ 2 {k\ ^ is required 
to ensure the traversability) . The properties of this wormhole are discussed by Dadhich et. al. 
[13]. 

It was shown by M0ller [16, 15] that the tetrad description of the gravitational field allows a 
more satisfactory treatment of the energy-momentum complex than does general relativity. We 
have then applied the superpotential method [16, 15] to calculate the energy of the gravitating 
system of the two solutions (21) and (24). As for the first solution we obtain E = m and there is 
no effect of the two parameters k\ and k 2 characterize the wormhole. As for the second solution 
the energy is E — 2m and this is due to the fact that the time-space components of the tetrad 

fields eo", e a ° go to zero as —= at infinity [35, 36]. The disappearance of the two parameters k\ 

V r 

and k 2 from the two Eqs. (41) and (43) may give an impression that these two parameters are 
not a physical quantities like the gravitational mass m. 



References 



[1] L. Flamm, Phys. Z. 17 (1916), 48. 

[2] M. S. Morris and K.S. Thome, Am. J. Phys. 56 (1988), 395. 

[3] M. S. Morris, K.S. Thorne and U. Yurtsever, Phys. Rev. Lett. 61 (1988), 1446. 

[4] J.D. Novikov, Sov. Phys. JETP 68 (1989), 439. 

[5] M. Visser, Lorentzian Wormholes-From Einstein to Hawking (American Institute of Physics, 
New York, 1996). 

[6] P. K. F. Kuhfittig, Phys. Rev. D,67 (2003), 064015. 

[7] T. A. Roman, Phys. Rev. D,47 (1993), 1370. 

[8] H. Epstein, E. Glaser and A. Yaffe, Nuovo Cimento 36 (1965), 2296. 

[9] M. Visser and D. Hochberg, "Geometric Wormhole Throats''' ',The internal structure of black 
holes and spacetime singularities, (Haifa, Israel, June- July 1997), gr-qc/9710001. 

[10] D. Hochberg, Phys. Lett. 5,251 (1990), 349; B. Bhawal and S. Kar, Phys. Rev. D,46 (1992), 
2464. 

[11] S. Kar, Phys. Rev. D,49 (1994), 862; S. Kar and D. Sahdev, Phys. Rev. D,53 (1996), 722. 

[12] D. Hochberg and M. Visser, Phys. Lett.,81 (1998), 764; Phys. Rev. D,58 (1998), 044021. 

[13] N. Dadhich, S. Kar, S. Mukherjee and M. Visser, Phys. Rev. D65 (2002), 064004. 

[14] D. Hochberg and M. Visser, " General dynamic wormholes and violation of the null energy 
condition", Advanced School on Cosmology and Particle Physics , Spain 22-28 June (1998), 
gr-qc/9901020. 

[15] C. M0ller, Mat. Fys. Medd. Dan. Vid. Selsk. 39 (1978), 13. 
[16] C. M0ller, Mat. Fys. Medd. Dan. Vid. Selsk. 1 (1961), 10. 
[17] D. Saez, Phys. Rev. D27 (1983), 2839. 
[18] H. Meyer, Gen. Rel. Grav. 14 (1982), 531. 

[19] K. Hayashi and T. Shirafuji, Prog. Theor. Phys. 64, 866, 883, 1435, 2222; 65 (1980), 525 

[20] F.W. Hehl, J. Nitsch and P. von der Heyde, in General Relativity and Gravitation, A. Held, 
ed. (Plenum Press, New York) (1980). 

[21] C. Pellegrini and J. Plebanski, Mat. Fys. Scr. Dan. Vid. Selsk. 2 (1963), no.3 

[22] K. Hayashi and T. Nakano, Prog. Theor. Phys. 38 (1967), 491. 

[23] K. Hayashi and T. Shirafuji, Phys. Rev. D19 (1979), 3524. 

[24] W. Kopzynski, J. Phys. A15 (1982), 493. 

[25] J.M. Nester, Class. Quantum Grav. 5 (1988), 1003. 



[27] V.C. de Andrade and J.G. Pereira, Phys. Rev. D56 (1997),4689. 

[28] V.C. de Andrade, L.C.T Guillen and J.G. Pereira, Phys. Rev. Lett. 84 (2000), 4533; Phys. 
Rev. D64 (2001), 027502. 

[29] H.P Robertson, Ann. of Math. (Princeton) 33 (1932), 496. 

[30] J.P.S. Lemos, F.S.N. Lobo and S.Q de Oliveira, Phys. Rev. D68 (2003), 064004. 

[31] N. R. Khusnutdinov, Phys. Rev. D67 (2003), 124020. 

[32] M.I. Wanas, Int. J. Theor. Phys. 24 (1985), 639. 

[33] C. M0ller, Ann. of Phys. 4 (1958), 347; 12 (1961), 118. 

[34] F.I. Mikhail, M.I. Wanas, A. Hindawi and E.I. Lashin, Int. J. of Theor. Phys. 32(1993), 1627. 
[35] T. Shirafuji, G.G.L. Nashed and K. Hayashi, Prog. Theor. Phys. 95 (1996), 665. 
[36] G.G.L. Nashed, Nouvo cimento B117 (2002), 521. 



